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ABSTRACT 


The in-line and lateral forces acting on cylinders and 
the in-line force acting on spheres placed in an harmonically 
oscillating flow have been measured. The periodic flow has 
meen Senerated in a2 U-shaped channel through the use of a 
pneumatic oscillator. 

Diewinaowmeiiereda and the Jift coefficients as well as 
the total force coefficient have been determined through the 
use Of Fourier analysis and found to depend on a period pa- 
mame ter . 

The results have shown that the transverse force acting 
on a cylinder in a periodic tfiow with zero mean velocity is 
as large as the in-line force. The frequency of the alter- 
nating transverse force depends on a period parameter and 
several frequencies may occur during a given period of oscil- 
Jation of the flow. 

It is recommended that the experiments be extended to 
Reynolds numbers as large as 10° in a similar but larger 


apparatus. 
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ieee ERODUCT ION 


A. EXPERIMENTAL JUSTIFICATION 

The subject of forces acting on bluff bodies immersed in 
time-dependent flows has been and will continue to be of in- 
meresct to fluid dynamicists, aerodynamicists, and practicing 
engineers for the special” reason that in nature neither the 
body nor the fluid which surrounds the body is ever ina 
state of steady motion. Even the flow behind a bluff body 
moving steadily through a fluid is accompanied by large scale 
unsteadiness. Thus, any type of unsteadiness of the ambient 
mow and/or the motion of the body introduces additional 
changes in the characteristics of the flow and its analysis. 

The drag and inertial forces are interdependent as well 
as time dependent, and the resistance coefficients obtained 
in unseparated flows are not applicable to Separated cases. 
Although indirect, the role of viscosity is paramount in that 
its consequences are Separation, vortex formation and shed- 
meee anda resultant alternations in the virtual mass. It is 
thus clear that it is necessary to determine the relation- 
ships between various resistance components in terms of the 
unsteadiness of the ambient flow and/or body motion, geometry 
of the body, the degree of upstream turbulence, past history 
me tne flow, etc. [Ref. 1]. 

However, the understanding of the behavior of time-depen- 
dent flows about bluff bodies may be advanced only by consid- 


ering the relatively more manageable cases and gradually 
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meeorating the information so obtained. This led to the 
present experimental investigation of the harmonically oscil- 


lating flow about circular cylinders and spheres. 


IS. SURVEY OF THE PREVIOUS INVESTIGATIONS 

A review of some of the previous investigations on time- 
dependent flows may best be presented by dividing the motions 
under consideration into various but admittedly arbitrary 
categories. 


1. Theoretical Analysis of Unseparated Time-Dependent 
Laminar Flows 


The study of the theory of separation-free time-de- 
pendent laminar flows has enjoyed particular attention, 
Pemely due to its practical significance and partly due to 
its relative mathematical simplicity. Surveys by Stewartson 
eer 2), otUart (Ref. 3S], and Rott [Ref. 4] reflect the cur- 
rent level of understanding. 

The unseparated class of unsteady flows, most of 
which results from the unidirectional or Ree ae accelera- 
tion of bodies in an infinite or bounded fluid medium, gives 
rise to an induced mass which must be added to the real mass 
of the body. In general, the virtual inertia tensor depends 
on the shape of the body, the nature of the fluid mediun, 
the orientation of the moving object with respect to the di- 
merron OL MOLIOn, and the depth of relative submergence 
from a free surface and/or solid boundary. 

The effects of viscosity and separation on added 


mass or the dependence of added mass of a given body on the 
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characteristics of the, time-dependent separated flow cannot 

yet be theoretically evaluated nor can it be experimentally 

separated from the total resistance experienced by the body. 
eae Ue libuae Vibration of Bodies in a Liquid 

Cuteuitoemotesesteand the Initial instants of Accel- 

Sue Onmermmoadtos trom a State of Uniform Velocity 

Vibratory motion has often been used [Refs. 5-8] as 
one of the experimental means to determine the virtual mass 
Peetficients. The technique has been restricted to ampli- 
tudes of motion which are so small that either separation 
does not occur or the characteristics of separation do not 
change due to the imposed, high-frequency vibration. Thus, 
the results obtained are not applicable to occurrences in 
which the duration of unsteady flow in one direction is long 
enough for separation EOmwoc Cur’. 

It is not necessary that the body be initially at 
rest to determine the added mass provided that one starts 
with a steady ambient flow about the body or with the body 
tone 2b constant velocity in a fluid otherwise at rest. 
Hamilton and Lindell [Ref. 9] have shown that the added-mass 
coefficient of a sphere determined by imposing a uniform 
acceleration onto its steady state motion is equal to its 
theoretically determined value as well as to its experimen- 
tally determined value obtained by small-amplitude, high- 
frequency oscillations. 

Rayleigh [Ref. 10] has oun that if the accelera- 
tion continues, io — V+ ae both the drag and the 


added-mass coefficients change with time. 
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The foregoing arguments suggest that the force ex- 
erted by the fluid on a body immersed in a time-dependent 
flow might be expressed as 


i BC, pA|V|V + C dV/dt + history-dependent i) 


Mf 


drag and inertia forces 


moewhtch A is the projected area of the body and M - the mass 


more the fPluid-displaced by the body. The drag coefficient Ch 
is to be taken equal to the quasi-steady state resistance 

coefficient and Cy to be taken equal to that measured or cal- 
culated from the potential theory when the value of the his- 


tory-dependent force is zero [Ref. 9]. 


3. Bodies Subjected to Large Amplitude Harmonic 
Oscillations ina Fluid Otherwise at Rest 


A number of significant studies |[Refs. 1i-i6] were 
made to determine the forces acting on bluff bodies under- 
Peine harmonic, large-amplitude oscillations in a liquid 
otherwise at rest. 

Odar and Hamilton [Ref. 11] measured the force act- 
ing on a sphere subjected to harmonic oscillations and pro- 
posed a force equation comprised of three parts: a steady 
drag component where the drag coefficient is the well-estab- 
lished coefficient for steady translation of a sphere; an 
inertial component whose coefficient is determined experi- 
mentally as a function of Wey arc Ca sac@eeleration; dad: diameter 
of sphere); and a history-dependent component whose coeffi- 


cient also is a function of V?/ad. 
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4. Bodies Subjected to Wave Motion with Zero Mean 
Velocity 


Considerable work [Refs. 17-23] has been done on this 
type of time-dependent flow in an attempt to predict the wave 
ces on piles and other submerged structures. In general, 
the total force acting on the body is assumed to be composed 
of a velocity-dependent and an acceleration-dependent force. 
“The resulting equation, ‘known as the Morison equation [Ref. 


16], is given by 


F = £C)pA|V|V + Cm, aV/dt (2) 


is 


where Oi = 1+C, C being the added mass coefficient. Experi- 
ments show that Ch and Ch depend on time and show consider- 
able scatter about their assumed mean values [Refs. 18, 21] 
over a cycle. 
Teme cales oubyected to Unidirectional Acceleration in a 

Fluid Otherwise at Rest or Unidirectional Unsteady 

Flow About Bodies Held at Rest 

substantial effort has been made to determine the 
Semponents of force acting on a body accelerating unidirec- 
Menally in wird ouherwise at rest. Iverson and Balent 
[Ref. 24] and Keim [Ref. 25] towed spheres, disks, and cylin- 
ders in unsteady motion through still water. It was con- 
@budaed that the total instantaneous force on these objects 
could be described by a single coefficient of the form 
C = oie AV Miietneae this, Coefficient was a function only Ct. 


a so-called acceleration modulus, ad/V*. Keim, however, de- 


tected an apparent Reynolds number effect for cylinders. 
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Laird et. al. [Ref. 26] measured the forces acting 
on accelerating cylinders and found strong evidence of devia- 
tion of the drag coefficient from the accepted values for 
Mmrcorm motion. They further found that the acceleration 
modulus did not correlate the resistance coefficient near 
boundary-layer transition. 

Sarpkaya and Garrison [Ref. 27] working with unidi- 
rectional flow with constant acceleration about cylinders 


and plates have found that both Ch and Ci depend on the his- 


h 
wemy OL MOtion, that the Morison equation (Eq. 2) is valid 
only for flows with constant acceleration, and that the 

acceleration modulus can correlate the data only if the ac- 


celeration is kept constant. 


Genome ecud Or Delfi—-Induced Transverse Oscillations of a 
peewee ntoa Finid ian Steady Motion 


Experiments have been directed to the effects that 
attend the oscillations of a body in the plane of the lift 
force. This type of interaction becomes particularly impor- 
tant when the oscillation of the structure is the necessary 
condition for the generation of the exciting forces. The 
most pertinent information from these experiments is that 
m@enlitt and draw forces act on the oscillating cylinder at 
the vortex shedding frequency and twice the shedding fre- 
quency, respectively, provided that the driving frequency on 
the cylinder is appreciably different from the shedding fre- 
quency. When the forcing frequency of the cylinder approaches 


the shedding frequency, the natural shedding frequency 1S 
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Mest and it ‘locks-in" to the forcing frequency. This syn- 
chronization persists over a range of frequencies which may 
be termed the "range of synchronization" [Ref. 28]. 
fee Unidirectional Oscillatory Flow About Bodies at Rest 

and obtveady Unatorm Viow About Bodies Subjected to 

wibeamwese Oscillations 

Relatively few studies have been conducted in this 
category. Chan and Ballengee [Ref. 29] examined the vortex 
emecaing from circular cylinders in an oscillating free 
stream. Their results suggested that "in an oscillatory 
freestream of 3 HZ and Reynolds number up to 40,000, the 


vortex shedding from a circular cylinder responds instantan- 


eously to the freestream variations" and that "the instan- 


taneous Strouhal number stays sensibly constant at 0.2 + 0.01." 


Hatfield and Morkovin |Ref. 30] studied the effect of 
an oscillating freestream on the unsteady pressure on a cir- 
@ular Cylinder. They have found that there is no significant 
coupling between the small-amplitude freestream oscillations 
and the vortex shedding. Their results would suggest that 
the drag coefficient associated with the mean flow would es- 
Seaeralily remain constant at its steady state value. On the 
other hand, Mercier [Ref. 31] who subjected cylinders to 
large streamwise oscillations found that the average drag 
coefficient significantly increases with fd/V and that the 
mite Of increase depends on the amplitude to diameter ratio. 
Combined, these results would suggest that the degree of 
coupling between the ineoqMencyecot oscillation and the vortex-— 
shedding frequency and hence the forces acting on the body 


strongly depend on the amplitude of oscillation. 
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Davenportcethet. S52} subjected bluff bodies to smail 
Pipettude Oscillations in a water flume and evaluated the 
drag and inertia coefficients through the measurement of the 
rate of damping of the amplitude of oscillations of the bod- 
1es, i.e. without measuring the forces acting on the bodies. 
Using the frequency parameter fa/V to correlate the data, he 


found large variations in Ch and Chr 


Se oCOPE OF THE PRESENT STUDY 

The present investigation is limited to an experimental 
mavestigation of the drag, inertia, and lift forces acting 
meme idrcular cylinders and the drag and inertia forces acting 
on spheres immersed in a periodically oscillating fluid with 


pemeoremean velocity. The fluid motion is characterized by 


V= Ne cos t C32) 


a 
T 
where ys 1s the maximum velocity, T the period of the oscil- 
Meoron, and V the instantaneous velocity. 

AS noted earlier, this type of fluid motion about cylin- 
ders and plates has been studied by Keulegan and Carpenter 
fRef. 18] through the use of standing waves in a rectangular 
Poin. ithe characteristics of the ambient fluid motion have 
been determined through the use of the appropriate standing- 
wave equations and the measured amplitude of the waves. 
meulesan and Carpenter determined the drag and inertia coef- 
merents and correlated them Rend cette so-called period pa- 
rameter Vala: They have found that the Fourier-averaged 


Values of the inertia and drag coefficients over a wave 
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cycle show considerable variations with Visa. HOR Te yh uncer s 
Vi /d equalling 15 was found to be a critical condition 
yielding the lowest value of the inertia coefficient and the 
mareest value Of the drag coefficient. For the plates the 
higher values of the drag coefficient are associated with 
the smaller values of Visa and the higher values of the in- 
ertia coefficient with the larger values of Vai /a- It was 
suggested by Keulegan and Carpenter that the parameter Vala, 
or equivalently, the relative fluid displacement A/d was the 
parameter of primary importance and the effect of Reynolds 
number was assumed to be of little importance on the values 
of the drag and inertia coefficients. It is noted, however, 
that the test method used by Keulegan and Carpenter allowed 
no contro] over the Reynolds number so that they had no con- 
venient method of testing this influence. Finally, it should 
also be noted that Keulegan and Carpenter did not measure 
wme litt forces acting on cylinders. 

Driscoll [Ref. 16] oscillated circular.cylinders in sim- 
ple harmonic motion in water otherwise at rest and found 
miu the Inertia coefficient 1s almost independent of the 
Reynolds number and highly dependent upon A/d and that the 
medrecOeclilicient 1S strongly affected by the Reynolds number 
and generally it decreases with increasing Re helds titi oes 

al 


Dreascoll's j{Ref. 16] Cy and Ch vVoluess are In error 


Since the length of the test cylinder was inadvertantly 
taken aS unity in the calculations instead of its actual 
Meneth of 1.375 inches. 
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The lift force averaged over the length of a circular 
pile was measured by Bidde [Ref. 33]. He found that the max- 


imum values of the lift coefficient expressed by 


= i 2 
Pimax 20%, oe 


LMAX So 
meaenes 60%, Of the longitudinal force for the rigid pile 
used and that Cr AX has a maximum average value Gm Jap OUlu. 1, On 
His data showed considerable scatter over the range of ei 
foe -eweneountercad in the experiments. Nevertheless, Cr MAX 
showed a better correlation with V_t/d than with the Reynolds 
number. This study does not shed much light on the under- 
Eeanaine® of the variation of the lift coefficient with either 
parameter since the lift varies along the pile in accordance 
fei the characteristics of the wave motion. 

Rance [Ref. 34] investigated the drag, inertia and lift 
forces acting on cylinders placed in an oscillating-flow 
water tunnel and found that the lift coefficient may reach 
values as high as 2.0. He did not attempt to correlate his 
results with the period parameter. Rance concluded that the 
lift force appeared to be significant compared with the sum 
P—piicmcrcaie andeinertia forces and that a combination of all 
the forces should be used in a practical design. 

Relatively few studies have been carried out on wave 
wiGeeseOn submerged spheres. Grace and Casciano [Ref. 35] 
measured wave forces on a subsurface sphere and was able to 


Solem ouly the drac coetficient which ranged from 0.3 to 


foot needats ana the method of calculation did not allow 
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them to obtain either a systematic variation of Ch With the 


period parameter or to evaluate Cure 
It is evident from the foregoing summary of the previous 
studies on unsteady flows in general and of the simple har- 
meme MOtions about cylinders in particular that the questions 
regarding the determination of the lift, drag, and inertial 
forces for various types of bluff bodies immersed in periodic 
flows remain largely unresolved. There are several reasons 
mem the difficulties encountered. Theoretically, the prob- 
lem is not manageable because of the unsteadiness of the 
flow and the indeterminate nature of the Separation points. 
Experimentally the difficulties stem from various sources in 
experiments. With oscillating bodies, the force resulting 
from the acceleration of the body must be separated by cal- 
culations. In the case of body immersed in an oscillating 
fmuid, the characteristics of the fluid motion have often 
been evaluated through the use of approximate wave equations 
rather than measuring them directly. Furthermore, either in 
mie Case of body oscillations or in the case of fluid oscil- 
lations, the motion was not entirely harmonic. This, as 
well as the ever present ocean currents, may have led to 
large scatter in the data presented in the references cited. 
Thus the values of Ci and Cy so derived are associated 
uniquely with the particular wave theory used to determine 
their numerical values. In general, attempts to use these 


values with other wave theories may result in large errors. 


famatiy, tt Should be mentioned that the values of Cys and Ch 
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which have been derived for cylinders placed horizontally in 
an oscillating fluid has often been identified in engineer- 
ing usage with similar coefficients obtained from tests re- 
lating to vertical cylinders. Obviously, the characteristics 
of the wave motion as well as the bottom- and free-surface 
effects will significantly affect the variation of these co- 
efficients with the significant parameters of the phenomenon. 
The foregoing reasons, along with related aspects of the 
problem gave rise to questions concerning the validity of 
the methods for the prediction of forces acting on bluff 
bodies immersed in periodic flows. The present investiga- 
tion is undertaken for the purpose of clarifying some of 
these problems through the use of a truly harmonic motion 
and to determine the drag, inertia, and the lift forces act- 
meon cylinders and the drag, and inertia forces acting on 


spheres. 
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Il. METHOD OF ANALYSIS 


The total force acting on cylinders and spheres in the 
direction of the ambient fluid motion is assumed to be com- 
prised of a drag force and an inertial force in the manner 
Similar to that first suggested by Morison et. al. [Ref. 17]. 


“For a cylinder it may be written as 


-_- i md? av 
Fo = 3C)pdg|Vv|V + Cy a7 Pf ae Go) 


mivere F is the force for the length 2, d the diameter of the 
body, oe the density of fluid, and V is the instantaneous 


velocity of the fluid. The Ch and C,, represent respectively 


M 


the drag and inertia coefficients. 
The corresponding equation for the sphere may be written 


as 


7a? dv 
Me noe at (6) 








2 
he C 


—|— i 
F aC p/P a 


in which F is the total force acting on the sphere. 
The velocity of the ambient flow is represented by 
V= ae EOc Ot (7) 


where Na denotes the maximum velocity, T the period, and 
Go = 27/T. The total force acting on the cylinder per unit 


length is in general given by 


F = £(t,T,V_,d,P,V) (8) 





Grouping the variables on the basis of dimensional rea- 


soning and introducing a = 27t/T gives 
a ee Vea 
. Vd 7 CE ae ee ar, (9) 


where V »t/y is a Reynolds number and Val/d;, which can also 
be Mressed as 2m7A/d, will be termed the relative displace- 
ment . 

Because of flow symmetry and the periodic nature of the 
force 

F(a) = -F(Catt) | (10) 

mis POSSible to express the force coefficient in a Fourier 
series. 


one a or Tere ole oO ft AcSin, ba + 


(11) 
TS neoOs Oa maCOS. Sast Beecos 5a + 


where the coefficients Ay. By are independent of @ and at 
most functions of Reynolds number and relative displacement. 


Fourier analysis may be used to determine the coefficients 


as: 
27 : 
ones ee en aN 
Ay 5 i oV2d da (22) 
O m 
SLigvek 
27 
Bp =lif ¥ocos Na ae (isp 
N TT 0 pvid 


Once obtained, the dependence of these coefficients on Rey- 
nolds number and relative displacement may be established 


provided the data are sufficient. 
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The general formulation, equation (11), may be reconciled 
with Morison's equation (5). Introducing V from equation 
(7) into equation (5) 


foe NG doe. noe 
ayzZad ~ a “yu am Sim @ - gr lcos ajcos o (14) 


m 


By the rule of Fourier 


211 
{ |cos a|cos a cos Nada 
a 0 
|cos alcos a = ar a 
f cos* Nada 
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where 
N+1 
er ae Zz, 8 
“oven oe “odd ou N(N -4)1 (79) 
The first three non-zero coefficients would then be: 
oy Svan. aa 2) OSI es ee 8/1051 (16) 
Introducing equation (16) into equation (11) with 
Ba = B,/ai 
B3 = B3-a3/a1(B,) 
Bee ey 2a 0b Gly) 
yields 
F _ 
pvid —Wiesil sat Assin 3a ~ AcSin 5a + 
m t J 
+ B,|cos ajcos a + B3}cos 3a]cos 3a 
fb | COS) 50) COS oa (18) 
Equations (18) and (11) may be compared. Writing 
v1 do Sin 3a Sei c 
a Se eee aa, 19 
4 @M Vv A; + As Sing ° “® sino 2 





Well 


i 1 ' 
Cy in B, Bacos of Bs cos 5a 
2 ae 


=a aS e e e 
[cos a] cosa [eos a]cos a — 


° t t 
Thus, if the coefficients A;, As and B3;, Bs vanish, the 


values of Ci and Cy remain constant for all phases of the 


cylinder motion 


Vat 
Ci = fac Aj G2) 


=n 


Substituting A, from equation (12) yields 


way ae a) I sin oda 
me wed |. yaa . 
0 ov7d 
m 
Ch = -2B) C235) 
substituting from equations (13) and (17) yields 

271 ne 

or = -AS F sin ada (24) 
0 pv d 


If the coefficients do vary with the phase a, the values 
piven by equations (22) and (24) are weighted averages. With 


this possibility in mind it is preferable to adopt 


t 
ad eosin & + Bi cos alcos a| + error (25) 
m 
or 
F T do ._. D 
as = — ——- aoa + 
avr , Cy, v sin o 5 | cos alcos a error (26 ) 
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where Ai, Bi, Ci and Ch are constant and error has the value 


CwEer = Tae, sineoe +toAs Sin 5a 


t 1 
+ B; cos 30 + Bs cos 5a 


This error may be obtained by subtracting the computed value 
of A; Sin a and B, | cos a|cos a from the observed F/pV_d. 

sathe .foreroing formulation was developed by Keulegan and 
Carpenter [Ref. 18] and equations (22) and (24) were the 
basis for the computer data reduction programs. 


For this purpose equations (22) and (24) were written as 


tollows 
eT ont i 
and 
me or Dit t 
Ch = SpdnhAz Ee COs ae) 6 Ga) (28) 


The corresponding coefficients for the sphere may be ob- 


tained through the use of equation (6) and through the use 





of a similar Fourier analysis. The results are given by 
pe eat ; Dig t 
and 
eo Qnt t 
Ch — Opn-a2Az eos T ) Ga) (30) 
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The error in force predictions or the difference betwecn 
the measured and the calculated forces acting on cylinders 
and spheres have been calculated not through the use of equa- 


tion (26) but rather through the use of the following expres- 


Sion: 
an SS - al 
Percent Error = ene ame e3t) 
mMeS“ max 
maewhich F i is the measured force and F is the calcu- 
mes Canlk 


Hated force. It is thought that the percent error calculated 
in this manner will be physically more meaningful than the 
predicted error given by equation (26). 

Computer programs have been devised for the evaluation 
Samune coefficients given by equations (27), (€28), (29), and 
(30). These programs are presented in Appendix A. 

Finally the maximum lift coefficient for each Vi /4a has 


been calculated through the use of 


FEMAX 


Ciomax = Epdgv" (32) 


2 


fm which PeOMAX is the maximum lift force in a given cycle 
for a given Vt/d. The corresponding computer program for 
Cr MAX is also given in Appendix A. No attempt was made to 
separately analyze the amplitudes of the harmonics of the 


iert force. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURE 


A. EQUIPMENT 

The basic oscillating flow system consisted of ‘ U-shaped 
vertical water channel. The cross-section of the test sec- 
tion was 18 by 20 inches. A photograph of the equipment and 
a schematic drawing of the channel are shown ae eicures 1 
and 2. The test bodies were mounted 2.5 feet below the free 
Surface of the still water level. The fluid in the U-channel 
was oscillated pneumatically through the use of a slider- 
crank mechanism which has periodically opened and closed the 
air supply line and also a large exit hole at the top of one 
Of the legs of the channel. A photograph of the oscillating 
mechanism is shown in fig. 3. 

The frequency of oscillation of the slider-crank mechan- 
ism was matched to the natural frequency of oscillation of 
the water column in the channel by means of reduction gears 
and through the use of a variable speed DC motor. The maxi- 
mum half-amplitude of oscillation of water column was 11 
inches and the natural damping of the oscillations was in 
the order of 1/8 inch per cycle without the use of the pneu- 
matic pulsating system. The water level at its minimum 
height in the test section was about 20 inches above the 
test body. The cylinders were manufactured out of plexiglass 
tubes or rods at desired diameters and at lengths approxi- 
mately 1/16 inch under the width of the test section. Self- 


aligning bearings were imbedded at each end of the cylinders. 
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Figure 2. Schematic drawing of the U-channel 
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The lateral and in-line force measuring devices consisted 
meevarious cantilever beams mounted to the outer surfaces of 
the test section. Eight piezoresistive strain gages were 
mounted on each cantilever beam and properly waterproofed. 
The drawing of one of the many Such force transducers is 
shown in figure 4. The end of each beam firmly fitted into 
mveeceli aligning bearings placed at each end of the cylin- 
der. | 

These transducers were repeatedly calibrated by hanging 
mera: at the mid section of the cylinders in the vertical 
@moenerizontal GQirections. Some of the force transducers 
were comprised of two pieces of thin cantilever beams cut 
orthogonal to each other. Thus they were capable of simul- 
taneously measuring both the lateral and in-line forces. 
Some transducers were capable of measuring only the lateral 
Sreethe in-line force. The one shown in figure 4 is of that 
by pe . It should be noted, however, that this particular 
transducer was rotatable + 90 degrees. This enabled one To 
measure first the in-line force and then, after a 90 degree 
mepation, the lateral force. In all cases both the in-line 
and lateral forces were measured at both ends of the cylin- 
fiereane compared with each other. In no case, did force 
mmves Geviate from each other, indicating a fairly uniform 
mesponse along the cylinder as far as the resultant forces 
are concerned. The size of the cylinders varied from 1.0 
inch to 2.5 inches. The spheres which were manufactured out 


Siealuminum, or hard plastics with diameters ranging from 
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feo tO 3.975 inches were hung with thin fish line from the 
force transducer placed at the top of the channel at the 
center of the test section. Attempts to prevent the spheres 
from swinging have resulted in the creation of undesirable 
oscillations in the force curves. Consequently, the spheres 
were allowed to hang freely and to respond to lateral forces. 
The amplitude of the swing naturally depended on the weight 
of the sphere. The maximum amplitude was approximately 0.25 
inch. As noted earlier, the fish line was connected to an- 
memer Lorce transducer, similar to the one shown in figure 
peeepeacea at the top of the channel. 

Throughout the investigation the monitoring of the char- 
acteristics of the oscillations in the U-channel was of 
peaie importance. As discussed in the introduction, most of 
m@emaiiticulties in the past in the determination of Cy Ch» 
and Cy MAX Cesuited trom the difficulty of creating a purely 
harmonic motion or from determining indirectly the character- 
forces Of the oscillatory motion. Even though the U-channel 
has provided, by its very nature, a perfectly Sinusoidal os- 
Seeerattven, both the instantaneous displacement and the accel- 
eration were continuously monitored. The instantaneous 
elevation in one leg of the channel was determined through 
the use of a capacitance wire connected to an amplifier- 
recorder system. Such wires have been used in the past to 
measure wave heights in open channels. The response of the 


merenvas found, through calibrations, to be perfectly linear 


Within the range of oscillations encountered. 





The instantaneous acceleration was measured by means of 
a differential-pressure transducer connected to two pressure 
taps placed vertically two feet apart at the mid-section of 
one of the walls of the U-channel. The instantaneous accel- 


eration was then calculated from 
AP = pra, C33)) 


where AP is the differential pressure, p the Alemenen On, ot id: 
& the distance between the pressure taps, and a, the instan- 
taneous acceleration of the fluid. The effect of the pres- 
sure drop due to the viscous forces over the length & was 
found to be Feces 

Evidently, the displacement and acceleration traces are 
in phase and may be used independently to calculate the ve- 
locity and the.displacement of the fluid. In fact such cal- 
culations have shown that the velocities and displacements 
calculated through the use of either traces did not differ 
more than 4%. Finally, it was gratifying to note that both 
the acceleration and displacement traces were nearly perfect 
Sine curves. Figures 5, 6 and 7 show the representative dis- 
placement, acceleration, in-line force and lateral force 
traces for the cylinders and figures 8 and 9, the displace- 
ment, acceleration, and in-line force traces for spheres. 
These traces which were obtained without any damping of the 
disturbances show the smoothness of the variation of the var- 
ious quantities and the degree of success achieved in obtain- 
me a purely harmonic motion. tees e= aso noted, from a 


brief perusal of the displacement and acceleration traces 
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displacement traces 


Force, 


Eremre 5. 








HE FO 


i 


ne 
i 
—. 

homme 
PN 
h 

‘ es 
; aH 


' 
ee eg tee —— 
, 


ae 
ae 


| 





t 
ee ee 
t 


displacement traces. 


Force, 


Pace 6. 


37 





' ; ' ‘ . 
' ace ' . a er Wh mah! | 
Qac ‘ ee ; : | ‘ 
Up eae hats 
i tees fp : i 
: eeedade Oar eens iste 


ons a 
fisbibaeedtscd: 


{ { fs 





i 





displacement traces. 


Bieure ~¢, Force. 
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displacement traces. 


Force, 


Figure 8. 
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displacement traces. 


Force, 


Figure 9. 
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iis the Gamping of the oscillations over a cycle 1s imper- 
wemeibly small when the fluid is let to freely oscillate in 


the channel. 


B. PROCEDURE 

From the traces similar to those described above, the 
amplitude of oscillation and the maximum velocity were cal- 
culated. Then the in-line force from the force traces was 
read and punched on cards for every 2.5 mm (this corresponded 
to At/T = 0.039645 for the period of T = 2.86 seconds). Then 
the drag and inertia coefficients were calculated through 
the computer programs given in Appendix A. This part of the 
procedure was identical for both cylinders and spheres. 

2ne calculation of the maximum lift coefficient consis- 
ted of the determination of the maximum amplitude of "lift" 
Sieetaceral force for a given Ve bad through the use of equa- 
men (32). 

It was realized before the start of investigation that 
not only the magnitude but also the frequency of the oscilla- 
b2ons of the lateral force is of major importance as far as 
the design of bodies subjected to harmonic fluid motions is 
concerned. With this idea in mind, the frequency or frequen- 
cies occurring in lateral forces in each cycle were deter- 
mined for each Vt /d. 

The maximum lift force in each cycle was also compared 
with the maximum value of the in-line force, as will be dis- 


cussed later. 
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Finally the phase angle between the occurrence of the 
maximum force and the maximum velocity was determined from 


the recorder traces for both the cylinder and sphere. 





IV. DISCUSSION OF RESULTS 


A. CYLINDER DATA 

The primary purpose of this investigation was the deter- 
mination of the forces exerted on cylinders and spheres in 
harmonically oscillating fluids. The secondary purposes 
were the@re-examination Of “the data obtained by Keulegan and 
Carpenter [{Ref. 18] and, whenever possible, the evaluation 
of the secondary effects such as non-harmonic oscillations, 
viscous forces, etc. on the forces obtained with harmonic 
oscillations. 

The drag and inertia coefficients for cylinders are 
shown in figures 10 and 11 as a function of the period param- 
eter Vai /d- The data follow in general the same trend of 
that obtained by Keulegan and Carpenter for very small values 


to eiia tomapoue 2.0 2and:C 


of Vi /d> e is nearly zero. 


M 
For Via /a equal to apeut 11.0, C 


D 
D reaches its maximum value 
Smeabpout 2.1 and Crs reaches its minimum value of about 0.8. 
Considering the fact that Cr = 1.0+k where k is the added 
mass coefficient, it is apparent that for some values of 

Vial /a the added mass coefficient may become negative. A 
Similar result has also been obtained by Keulegan and Carpen- 


ter. For larger values of Vall d: Comers nearly equal to 1.2 


M 
and Ch approaches 1.4, within the range of Vzt/d values en- 
countered in the present investigation. 


Mean lines have been passed through the data shown in 


figures 12 and 13 and these mean values are compared in 
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figures 14 and 15 with the corresponding mean lines obtained 
by Keulegan and Carpenter. Evidently, for smaller values of 
Vint /d the results of the two independent studies are compar- 
able. For larger values of Vial /d: however, the Keulegan and 
Carpenter study shows relatively larger values for Cire Even 
though it is not easy to explain the difference, it may be 
conjectured that the relatively few data points obtained by 
Keulegan and Carpenter through the use of a 0.5 inch diameter 
cylinder may have been in error or the standing waves corres- 
ponding to the large values of V b/d may not have followed 
the particular analysis used in the prediction of their char- 
Meweri’sStics. Furthermore, in this range of Vat/d values, the 
mertaal force is small relative to the drag force. This in 
turn leads to an unstable determinant in the equations yield- 
ing Cn and Ch: 

It has been .conjectured first by McNown and Keulegan 
meet. 21] and later by Sarpkaya and Garrison [Ref. 27] that 
there should be a unique relationship between Ci and Ch allele 
each of these coefficients are dependent only on Vt/d- oO 
explore this matter further, a plot has been prepared as 
shown in figure 16. The numbers on this figure give the 
corresponding V/d Villec aBVven unOuUc hon NO particular 
@ractical Significance, the variation of Cy with Ch is indic- 
ative of the existence of a complex relationship between the 
two coefficients. The data points are a step further re- 


moved from the actual fluid motion through the elimination 


of Vial /d: Thus, such a plot may be of some significance if 
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and when the problem could be analyzed numerically or other- 
wise. 

The phase angle between the occurrence of maximum in-line 
force and maximum velocity which is also a measure of the 
ratio of the relative magnitudes of the inertia and drag 
forces, is evaluated from the traces and is shown in figure 
17. Evidently for ae values less than about 8 or 10 
there 1S a large phase pe eens between the occurrence of 
the maximum force and the maximum velocity. In this range, 
the maximum velocity leads the maximum force. For Vn /d 
values larger than about 20, there is a very small phase dif- 
ference between the two quantities cited. For very large 
values of V/d:, the phase angle drops nearly to zero. It 
may thus be concluded that the inertial forces dominate the 
motion for Vial /a values less than about 8. For larger values 
of Viai/d; the drag force is dominant. Thus for structures 
Subjected to wave forces with Viyt/d values less than about 8, 
the drag forces may be ignored. The alternate is true for 
Vial /d larger than about 10. 

The maximum lift coefficient, defined by equation (32), 
is shown in figure 18 again as a function of Vai /d- Such 
data, which are presented herein for the first time, have 
several significant features. Firstly, the maximum lift co- 
efficient may reach a value as high as 3.0. Secondly the 
lift coefficient has several maxima and minima dependent on 
Vial /d- jie 1ift traces show that the lift force does not 


begin to develop for Vint /d less than about 4.0. Then the 
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Min cOeCcificient increases rapidly reaching its first maxi- 
mum value of about 2.85 at Vt/d =O 02e tne SigniitTeance 
Seine existence of large lift forces between 10 > Vt /4a > 4 
lies in the fact that even when the forces in line with the 
fluid motion are essentially of inertial nature (V t/a <—20) 8 
there is still a region, (4 < Vt/d < 8) in which consider- 
able feteral EObhececumay act wpen the cylinder. . Thus, if one 
were to ignore separation and vortex shedding and determine 
through the use of inviscid fluid flow analysis the forces 
acting on the cylinders, the results should not be extrapo- 
lated beyond Vani /a > 4, 

The lift coefficient sharply decreases in the range 
10 < Vt /d < 15 and then increases again to a value of about 
3.0. Subsequently, Cr MAX gradually decreases to about unity. 
It is apparent from the data shown in figure 18 that there 
i considerable scatter in the range of Vt/4 from about 18 
to 25. The reasons neither for this scatter nor for the rap- 
mararop and rise for Vial / a from 10 to 18 can be explained in 
Simple terms. They can, however, be attributed to the com- 
plex interaction of the shed vortices aft and fore of the 
cylinder and to the fact that we have not decomposed the lift 
force into its harmonics. 

Several attempts have been made to visualize the flow 
with various dyes and particles and to obtain motion pic- 
tures. These efforts were only partially successful in iden- 
tifying the instantaneous positions of the vortices. Efforts 


were also made to visualize the flow pattern in a water 
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table with a free surface. The strong effects of the surface 
tension have prevented the identification of the vortices. 

At small values of V fd beawtatul patterns of acoustic 
streaming have been observed. 

As cited in the development of the so-called Morison's 
equation, the coefficients Cu and Cy are Fourier averaged 
and thus assumed constant over a cycle. imeeealeaty, thie wis 
steady nature of the flow renders these coefficients time 
dependent and thus gives rise to some differences between the 
measured and calculated forces. This difference, called er- 
ror or remainder force, could be estimated either through 
the use of the additional terms appearing in equation (11) 
or through the evaluation of the percent difference between 
the measured and calculated forces. The former technique 
which was used by Keulegan and Carpenter [Ref. 18] may still 
differ from the actual error between the measured and calcu- 
lated forces. With ae view in mind it was preferred to 
calculate the actual error in each cycle for representative 
values of et [a as well as the error corresponding only to 
pace maximum force for all values of Vint /d- For this purpose 
the percent error was written as 


ie - F 


d <i a (34) 
MES «mas 
and 
vo Pines(max) 7 eee (35) 


Fnes(max) 


o6 





The values of A for representative values of Vai /d are pre- 
Sented in Appendix B. The variation of A* with V w/a ie 
Shown in figure 19. Evidently the maximum value of A* occurs 
at Vat/d =~ 12 and is in the order of 15%. Thus, the designer 
must multiply the calculated values with A* in determining 
the maximum actual force acting on the cylinders. 

So far attention has been paid to the calculation of the 
lift, drag, and inertia forces and nothing has been said 
either on the frequency of the lateral or lift forces or on 
the ratio of the maximum lift force to the maximum in-line 
lorce. 

As cited earlier, the frequencies of lateral forces in 
each cycle for each Vint /a have been evaluated for all cylin- 
ders tested and plotted in figure 20. Such a plot which has 
never been presented in the literature before, shows several 
interesting features. Firstly the vortex shedding frequency 
is not a constant fraction of the oscillation frequency 
(f = 1/T) and that several frequencies may occur during a 
given cycle. This may appear to be a priori evident consid- 
ering the fact that in an oscillating flow the instantaneous 
velocity of fluid varies from zero to ve and that there can- 
not be a single Strouhal frequency. Secondly, there does 
not seem to be any lateral force for Vi /a less than about 
4.0. Considering the fact that Vi /a may also be regarded 
as a measure of the relative displacement of a fluid particle, 
for small relative displacements vortices are shed at two 


Gistinct frequencies: one at the oscillating flow frequency, 
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memes, and the other at 2f. For some values of Vit/d, the 
vortices may be shed at as many as three or four different 
frequencies, namely, at fT aetnein Ste eeLe,, Kor Waly lar- 
ger than about 35 the vortex shedding frequency is limited 
to a single frequency of 2f. 

The significance of the existence of multiple frequencies 
mumperiodic flows stems from the fact that a structure may, 
under certain circumstances, come into resonance with the 
Mmeeguency of the exciting lateral forces acting on it. For 
waves of large periods this frequency may be sufficiently 
small to match the natural frequency of highly elastic thin 
members. It is not expected that highly stiff and large 
cylindrical structures would be excited by long waves. Thus 
aeevylindrical pile may go into resonant motion at relatively 
smaller wave periods at or near the frequencies twice the 
frequency of waves. The subject of vortex-synchronization wil] 
not be discussed here further Since the cylinders were held 
almost rigidly at rest in the pulsating flow. 

As to the magnitude of the lateral force relative to the 
maximum in-line force, a plot has been prepared representing 
this ratio in terms of Vi/d. The purpose of this effort 
was to emphasize the fact that the lateral forces could, un- 
der certain circumstances, be considerably larger than the 
in-line forces. In fact in the range 5 < Voto < 15, the 
Poaderatio is as large as 1.5. Thus in calculating the max- 
imum total force acting on the cylinder one must consider 


mnie svectorial sum of the lateral and in-line forces. With 
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this view in mind a total force coefficient defined by 


] 2 4 Fe 
P BD max ) FL (max) 


1 Zz 
zoU_ dk 


C 


’ (36) 


was calculated for two representative cylinders and the re- 
sults were plotted in figure 21. Evidently, Cr may reach 
Waves aS high as 4.0 in the vicinity of Vyzt/d a. A Oe ae 
Should be noted that the total force coefficient calculated 
in this manner is somewhat larger than the actual total force 
coefficient since the maximum lift force does not necessarily 
occur at the same instant as the maximum in-line force. How- 
ever, considering the complexity of the lift force and the 
possibility of the simultaneous occurrence of the maximums 
of the two forces it was preferred to calculate the total 
force coefficient as given above. No attempt was made to 
analyze in detail the phase angle between the first and high- 
er harmonics of the lift force and the maximum in-line force. 
A perusal of the data has shown that there are many instances 
where the two forces occur almost Simultaneously. It is evi- 
dent from the foregoing discussion that what is most impor- 
tant for the design of structures in such flows is not the 
drag and inertia coefficients for the in-line force but 
mmener a total force coeificient which fully accounts for 
both the in-line and the transverse forces. 

The dependence of or and Ch on the Reynolds number was 
mivesticated by plotting both coefficients as a function of 
Reynolds number for various values of Val /d.- As seen from 


mtioures 22 and 23 neither of the two coefficients depend 
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upon Reynolds number within the range of Reynolds numbers 


investigated. 


B. SPHERE DATA 

The inertia and drag coefficients calculated through the 
use of equations (29) and (30) are presented in figures 24 
mmemezo. Lhe inertia coefficient starts at its theoretical 
value of 1.5 and gradually drops to about ee hor Vat /d=15. 
The sphere data show considerably less scatter than those 
for the cylinder. Apparently, the added mass coefficient of 
a sphere in periodic flow could drop below its theoretical 
value of 0.5 for Ved larger than io. 

ime Grag coefficient shown in figure 25 gradually rises 
Pemabout 0.8 at Vi /d = 18 and then decreases very slowly 
with increasing values of the period parameter. Comparison 
of the Crs and Ch values with those obtained with the cylinder 
Shows that whereas in the case of cylinder Cr exhibits a min- 
imum at about Weel ad = 12, in the case of spheres, it de- 
creases gradually, as noted above. On the other hand, Cy 
values follow similar trends. 

The dependence of Crs and Ch on the Reynolds number was 
investigated by plotting both coefficients as a function of 
the Reynolds number for various values of Val/d. (See fig- 
ures 26 and 27). Apparently, neither of the two coefficients 
depend upon the Reynolds number within the range of Reynolds 
numbers investigated. 


The percent error A* given by equation (35) was calcu- 


lated for the sphere and is presented in figure 28. 
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Evidently the maximum value of A* may be as high as 15% at 
Vn /d Values of about 12 and 22. The reasons for the occur- 
rence of the maximum errors at the particular Vani /d values 
is not easy to explain. The first maximum takes place at 
the Vnt/d value where both the drag and inertial forces are 
quite significant and the maximum force is nearly 100 degrees 
out of phase with the maximum velocity as seen in figure 29. 
The second absolute maximum of A* occurs at a Vat/d value 
where the maximum force and the maximum velocity are nearly 
in phase. 

mmaliy, the variation of Ch with Cys 1s shown in figure 
30. The data show a trend similar to that presented for the 
cylinder. Although, these coefficients are averaged over a 
cycle and thus do not represent their instantaneous varia- 
Tions over a cycle, the existence of a unique relationship 


between them is evidenced by figure 30. 
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V. CONCLUSIONS 


The data presented herein warrant the following conclu- 
mions: 

1) A U-shaped channel is ideally suited for the genera- 
tion of perfectly harmonic flow oscillations about bluff 
bodies. Even though the period of the oscillations must be 
kept constant, the variability of the amplitude of oscilla- 
tions and the size of the test bodies enable one to vary the 
basic parameters of the investigation. Furthermore, some of 
the difficulties (e.g. undesirable vibrations, the separa- 
tion of the inertia of the test body, etc.) encountered in 
the experiments conducted by oscillating the test bodies are 
completely eliminated. 

2) The drag and the inertia coefficients obtained for 
the circular cylinder follow in essence those obtained by 
Keulegan and Carpenter. The present data are considered to 
be more accurate and reliable, particularly for the inertia 
coefficient, than those obtained by Keulegan and Carpenter. 

3) The finding that there is no noticeable correlation 
between the Reynolds number and the drag and inertia coeffi- 
cients confirms the similar conclusion reached by Keulegan 
and Carpenter within the range of Reynolds numbers encoun- 
tered. 

4) The above conclusions which mainly concern the con- 
firmation of Keulegan and Carpenter's historic data raise 


the important question that the scatter noted in the Cu and 
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Cy values obtained from the experiments conducted in ocean 
and noted by Wiegel [Ref. 36] may not be entirely due to the 
difficulty of the experimentation or due to the experimental 
errors. The scatter may rather be due to an additional fac- 
tor not present either in Keulegan and Carpenter's work or 
in the present work. It may be conjectured that the effect 
of the ever present ocean currents on wave-generated oscil- 
latory flow and possibly very large Reynolds numbers, for 
which the flow may become critical and the boundary layer 
about the body turbulent, may be responsible for the observed 
scatter. This aspect of the problem is currently being in- 
vestigated through the use of an oscillating flow water tun- 
mel. 

5) The most important and previously unexplored si mekom mares 
of the present investigation is that the lift or the trans- 
verse force is as large as the in-line force. The frequency 
of the alternating force depends on Vit/a and various fre- 
quencies may occur ina given cycle. These frequencies may 
occur at the frequency of the flow oscillations and at two, 
three or four times the flow oscillations. It is also pos- 
Sible that additional harmonics may be superimposed on these 
oscillations. These findings make it mandatory that the 
forces acting on cylinders be calculated by considering the 
vectorial sum of the in-line and transverse forces. One 
must also consider various DrCatenedes of oscillation in 
eeecane the fatigue characteristics and the structural in- 


tegrity of piles. 





oti tey Of eEne Sipnitaeant Contribution of the lift 
forces to the total force one must also raise the question 
as to whether undue amount of effort and attention have been 
given in the past to the determination of the components of 
the in-line force and to the use of the Morison equation for 
the prediction of the forces acting on cylinders when the 
transverse force iS as much as or larger than the in-line 
force. It appears from the discussion of the total force 
coefficient that there is a fairly good correlation between 
moe total force coefficient and Vini/d- It also appears that 
the largest value of the said coefficient occurs at V_t/d=10. 
The use of such a coefficient together with a knowledge of 
the frequency of lateral oscillations may be more meaningful 
than the use of the Morison equation which deals only with 
the prediction of the in-line forces. 

7) The force coefficients for the sphere follow in gen- 
eral the same trends as noted for those for the in-line 
force acting on cylinders. No attempt was made to measure 
the lateral forces acting on Spheres. Apparently, these are 
more random in nature and thus more difficult to couple with 
the in-line forces to determine a total force coefficient in 
terms of Vjt/d. The sphere data, like the cylinder data, 
show that there is no correlation between the force coeffi- 
cients and the Reynolds number within the range of Reynolds 


numbers investigated. 
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VI. RECOMMENDATIONS FOR FURTHER STUDIES 


In view of the success achieved with the oscillatory flow 
channel used in the present investigation it is recommended 
that the experiments by extended to higher Reynolds numbers 
(a0proximately 10°) in a similar but larger channel and the 
in-line and transverse forces’*be measured and analyzed in 
the same manner for cylinders, spheres and plates placed not 
only in the middle of the test section but also close to one 
of the walls of test section. The latter type of experiments 
will increase our understanding of the forces acting on 


large pipe lines, placed on or near the ocean bottom. 
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APPENDIX A 


Computer Programs 


ee eo CCC CCC CECCC CCE COCECC CECE CECE 
p AVERAGE CD AND CM CALCULATIGNS FGR CYLINOERS 2 
C TIME=CIMENSICNLES TIME (TIME/PERIOC) C 
C BETA=DIMENSIONLES DISPLACEMENT (CUMAX#PER/DIA) c 
C CLECYLINDER LENGTH IN FEET c 
C DIA=DIAMETER GF CYLINDER IN FEET g 
C AMP=AMPLITUDE OF MOTION IN FEET C 
C PER=PERIOD IN SECS (CHART PERIOD/CHARTSPEED) c 
C Z2eGORes COEFFICIENTS c 
C CM=INERTIA COEFFICIENT C 
C CuSonno COEFLICIENT C 
C REMF=REMAINGER FUNCTION cE 
C Witsoe =-F UURLTER COEFFICIENTS c 
c N=NUMBER GF DATA SETS C 
5 CN=CGNVERTIGN FACTGR ‘ 
munmmnewme mee CCeCCCCCCCCECCCCCCCCECCCCCCCCCCCCCCCCCCCCCCCCCCCC 

G=32.1/4 

RROC=6 2e4/G 

PT=3.14159 

CNU=0.0000105 
C CL=1.4946 
DNETIATE DATA SETS 

N=] 
P DC 100 [=1,N 
4 READ IN PARAMETERS WHICH CHARACTERIZE THE DATA SET 
: READ(5,10)01A,AMP,PERsNCARD,UMX,CN 
: CGneUuUTE FORCE COEFFICIENTS 

ZL=Ef2*PER**2/(PT**S3*DTA**2*CLYERHOXAMP ) 
: Z2=-3*P ER¥*2/(8*RHO*DIA*PI*XCL*AMP*X*2) 
CCMPUTE BETA 
c BETA=2*PI*AMP/DIA 
: CCMPUTE REYNOLDS NUMBER 
' REYNO=(UMX*DIA)/CNU 
c PRINT OUT PARAMETERS AND COEFFICIENTS 

WRITE (6,15) 014s AMPsPER,2Z12Z22,3UMX—CN 

WRITE(6,20) 

TIME=0.0 

CM=0.0 

CD=0.0 
Ps DELTAT=0.034965 
: START DATA REDUCTION 

00 290 J=1,NCARD 

READ(5;,25)F 

F=CN*F 

ALPHA=2*PI*TIME 

SINA=SIN(CALPHA) 

CuSA=COS( ALPHA) 

FSINA=DELTAT*F*SINA 

FCGSA=DELTAI*F*COSA 

CM=FSINAtCM 
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CCCCCCCCCCCC CCCCCCCCE CECE CECCCCECCCC EEC ECCCCC CECE CCCCCC CCC 
C AVERAGE CD AND CM CALCULATIONS FCR SPHERES C 
C TIME=DIMENSIONLES TIME (TIME/PERICD) C 
C BETA=DIMENSIONLES DISPLACEMENT (UMAX*®PER/DIA) C 
C DAO LAHETER GE SPHERE IN FEET é 
G AMP=AMPLITUDE OF MOTION IN FEET C 
C PER=PERTOD IN SECS (CHART PERIOD/CHARTSPEED) C 
C ZI=FORCE COEFFICIENTS E 
C CM=INERTIA COEFFICIENT ie 
C CE=-ORAG COEFeICTENR é 
C REME=REMAINDER FUNCTION C 
C Pi Pil-=—BUR TER COEEEICTENTS C 
C N=NUMBER OF DATA SETS C 
C CN=CGNVERTION FACTOR C 
Mme (CEC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
G=32.174 
RHG=62.4/6 
P1=3.14159 
: CNU=0.0000165 
C iQmhATe DATA SETS 
: N=1 
C READ IN PARAMETERS WHICH CHARACTERIZE THE DATA SET 
: DO 100 I=1;N 
: READ(5; 10) DIA yAMP,PERyNCARD,UMX,CN 
C CaMENTE FORCE COEFFICIENTS 
7 21=3*PER*#Z/ ((P1¥*3)* (OL A*3) XRHGXAMP ) 
, 7 2=—3%P ER¥¥2/ (2*QHO*(PI*2)*(DIA*2)*(AMP**2)) 
“ GEMPULE BETA 
: BETA=2*PI*AMP/DIA 
COMPUTE REYNOLDS NUMBER 
: REYNO=(UMX*D1A)/CNU 
C PRINT GUT PARAMETERS AND COEFFICIENTS 
WRITE(6é,15)D0DIAsAMP;PER;Z213Z2,UMX;CN 
WRITE(6,20) 
TIME=0.0 
CM=0.0 
CD=0.0 
; DELTAT=0.034965 
. START DATA REDUCTION 
DG 200 J=1,NCARD 
READ( 5,25) F 
F=CN*F 
ALPHA=2*P1¥*TIME 
S INA=S IN( ALPHA) 
COSA=COS{ ALPHA) 
FSTNA=DELTAT¥F#SINA 
FCOSA=DELTAT*F*CGSA 
CM=FSINAtCM 
; CD=FCOSAtCD 
C PRINT OUT COMPONENTS FOR DRAG AND INERTIA COEFFICIENT 
NRME C64) TIME sALPHAyGUSASS INA, F,ECOSAs FSINA 
TIME=TIME+DELTAT 
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